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Heat production and dissipation induced by current flow in nanostructures is of primary impor- 
tance to understand the stability of these systems. These effects have contributions from both 
electron-phonon and electron-electron interactions. Here, we consider the effect of the local electron 
and ionic heating on the conductance of nanoscale systems. Specifically we show that the non-linear 
dependence of the conductance on the external bias may be used to infer information about the 
local heating of both electrons and ions. We compare our results with available experimental data 
on transport in D2 and H2 molecules. The comparison between experiment and theory is reasonably 
good close to the lowest phonon mode of the molecule, especially for the D2 molecule. At higher 
biases we cannot rule out the presence of other effects like, e.g., current-induced forces that make 
' the scenario more complex. 
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INTRODUCTION 



The idea of building electronic devices from nanostructures has gathered a lot of attention due to the high expec- 
tations in terms of size reduction and power dissipation Encouraging progress has been made in experimental 
techniques and theoretical modeling towards this aim However, a fundamental and technologically important issue, 
5— 1 ■ namely local heat production and dissipation in these systems has attracted much less attention [1, 0; H B S E3] • 
It has been argued that since the electron inelastic mean free path is large compared to the dimensions of a 
nanostructure, no energy dissipation occurs in the nanostructure region. However, nanoscale systems carry very large 
current densities compared to bulk electrodes. This implies an increased number of scattering events per unit time 
and unit volume so that interactions among electrons or among electrons and phonons are particularly important. 
\ In addition, the reduced size means a small heat capacitance: any small energy transfer from the current-carrying 
electrons to local ionic vibrations or other electrons in the system may induce a substantial heating of the nanostructure 

o: 0. 

So far, direct measurements of the amount of energy locally dissipated in a nanoscale system have been beyond 
our reach. However, new experiments have considered the indirect effects of local heating on accessible quantities 
-H ■ @, @, EH, EH- For example, in 0, 0| an effective ionic temperature is determined via the force needed to break 
the chemical bonds between molecules and the adjacent leads. These experiments indirectly probe the local ionic 
temperature, the contributions due to electron-electron interactions, and corresponding local electron heating @, @|- 
Here, we discuss another possible indirect method to probe both the local ionic and electron temperatures via 
the non-linearities in the DC conductance of nanostructures. We will compare our results with the experimental 
conductance of simple molecules such as D2 and H2 sandwiched between two Pt leads as studied in Ref. [l3| (and 
■ references therein). 

In order to address the above issues we need a theory that takes into account both energy production and dissipation 
on an equal footing. A full quantum-mechanical description in terms of many-body states for the present non- 
equilibrium problem seems hopeless. Instead, we have previously shown that a much more "economical" hydrodynamic 
theory in terms of the single-particle density and current density may be derived for nanostructures In this paper, 
we first review such a theory, and later on use it to study the effect of heating on conductance. 



CLASSICAL HYDRODYNAMICS 



In the following, we will refer to some concepts of classical hydrodynamics. For completeness, we repeat here some 
of those concepts, while a more comprehensive description of the dynamics of classical fluids can be found in many 
textbooks 

The dynamics of a classical viscous fluid is usually described by the so-called Navier-Stokes equations for the 
single-particle density, n(r,t), and the velocity field, v(r,t), (ratio between the current density and the density) 



D t n(r,t) — — n(r, t)V • v(r, t) 
mn(r,t)D t Vi(r,t) = -ViP(r,t) + V^^r,*) - n{r,t)ViV ext {r,t) 



(1) 
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where P(r,t) is the pressure, 7tij(r, t) is the Navier-Stokes stress tensor 



VjVi(r,t) +ViVj(r,t) - -SijV ■ v(r,t) 



(2) 



and V ext (r, t) is the external potential. [Throughout the paper, V, corresponds to the derivative with respect to the 
7-th spatial component (i = {x, y, z}), and summation over repeated indexes is understood.] In these equations, the 
operator D t = d t + v(r,t) ■ V is the so-called "convective" derivative, while the viscosity coefficients, n and (, are 
the shear and bulk viscosity of the liquid, respectively. The viscosity coefficients have their origin in the approximate 
nature of the Navier-Stokes equations and in the particle-particle interaction 17 1. The first equation in |T]) is the 
continuity equation and states the mass conservation when sources or sinks are not present. The second equation in 
([1]) is the force equation: the left hand side is the acceleration of a small volume of liquid subject to the internal forces 
(due to pressure, particle-particle interactions) and external forces (F ex t — — W ex t)- It is important to realize the 
approximate nature of these equations: in classical physics the very basic concept of particle density n has a meaning 
only in a coarse-grained sense, i.e., with respect to volumes of the liquid small compared to the other relevant scales 
of the problem, but large to contain "enough" particles so that a continuum mechanics can be developed. In the 
opposite condition, one has to revert to the solution of the Newton equations of motion for each particle. Due to the 
continuous spatial nature of wave- functions, the above limitations do not pertain to Quantum Mechanics, for which 
a hydrodynamic description can be formulated exactly. 



HYDRODYNAMICAL FORMULATION OF QUANTUM MECHANICS 

Ever since the formulation of the Schrodinger equation of motion for complex wave- functions, there have been several 
attempts to formulate Quantum Mechanics in terms of classical quantities. The degree to which these attempts have 
been successful is still undecided, since the use of words like "particle" , "trajectories" , and "directions of propagation" 
is widespread in the modern scientific literature. One such attempt was made at the dawn of Quantum Mechanics, 



in 1926, by Madelung [17J, [18| who showed that the Schrodinger equation for a single particle is exactly equivalent to 
a set of equations of motion for the particle density and "velocity" . For this single-particle problem, the velocity is 
defined as the variation of the phase of the wave-function with position, and thus seems a mere mathematical tool 



171 ] . An equivalent, but more transparent definition of this velocity field is 



v(r,t) = ——-, (3) 
n(r, t) 

where j(r, t) is the current density. This definition is valid for the points r for which n(r, t) ^ 0. It is remarkable that 
the equation of motion for this velocity is governed by the external forces, plus a "quantum mechanical" contribution, 



known as "Bohm stress tensor", that has not a classical counterpart [T7[. Indeed, if we start from the Schrodinger 
equation for the wave-function, of a particle in the presence of the external potential V ex t, (h = e = 1 throughout 
this paper, where e is the electron charge) 

id t V(r, t) = _-!v 2 *(r, f) + V ext (r, i)*(r, t) (4) 

we can rewrite ^(r, t) in terms of two real functions of time and position, R and 5, as 

*(r,t) = R(r,t)e iS ™. (5) 

It is a simple exercise to show that, if one defines the density n(r, t), and the velocity v(r, t), 

n(r,t) = R 2 (r,t) = \*(r,t)\ 2 , (6) 
VS(r,t) 

v(r,t) = , (7) 

771 

then the equations of motion 

dtn(r,t) = -V- [n(r,t)v(r,t)], (8) 
md t v(r, t) = J-V ( V p ^ r '^ ) - «77j(r, t) • V7j(r, t) ~ VV ext (r, t) (9) 



2/77 V R{r,t) 
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hold. Eqs. {8]| and ([9]) have a clear physical interpretation: The quantum mechanical system is equivalent to a fluid 
whose dynamics is governed by the Euler equation ([9]) subjected to the force exerted by the external potential [ll|, 
and an internal force whose origin is purely quantum mechanical. 1 Moreover, the dynamics conserves the mass, 
i.e., the total probability, and then the continuity equation ([5]) holds [l9j| . The solution of the equations of motion 
([5]) and © is equivalent to the solution of the Schrodinger equation. It is interesting to point out that the quantum 
mechanical force can be expressed in terms of the Bohm stress tensor, 

Pi,j = ~^-n(r, t)ViVj \n(n(r, t)) (10) 
4m 



and 

r 

' 2m "V R{r, t) 

If one introduces the convective derivative the equations of motion ([H]) and ^ assume the well known form of the 
Navier-Stokes equations of motion 

D t n(r, t) = -n(r, i)V ■ v(r, t), (12) 
mn{r,t)D tVl {r,t) = -VjP id {r,t) - n(r,t)\7,V ext (r,t). (13) 

Eqs. (fi"2"|) and (TIB")) are formally identical to the Navier-Stokes equations (Tl]) for a classical fluid. However, unlike the 
Navier-Stokes equations which describe an approximate dynamics of the many-body classical fluid, Eqs. (| 1 2[) and (fl3|) 
are exactly equivalent to the Schrodinger equation: no approximation has been made in their derivation. 

While this approach to Quantum Mechanics may appear as a simple attempt to recover a classical description of 
quantum phenomena, over the years it has proven to be a very useful tool to describe the dynamics of quantum 
systems in several contexts ranging from condensed matter physics to nuclear physics (see, e.g., (20j . and references 
therein). More recently, we have shown that a hydrodynamic description of the electron flow in nanoscale systems 



leads to theprediction of novel phenomena, like the existence of a dynamical (viscous) resistance [21|, turbulence 
[3 0, HE 24 1 , and local electron heating and its effect on ionic heating @, 0] . 



Here, we describe our hydrodynamical approach to transport in nanostructures. As a first step, we need to generalize 
the derivation of the equations of motion (JT3J)- (U3J) to the case of a many-body interacting system. We follow closely 



the formalism presented in Refs. [lj, |25j. (See also Ref. [26j for a general formulation of the dynamics of a many- 
particle electron system.) We describe the dynamics of the system via a field creation (annihilation) operator rp'(r,t) 
(ip(r,t)) which evolves in time following the Heisenberg equation of motion 

id t i>{r,i) = --Lw 2 ^jj(r,t) + V ext (r,t)^(r,t) 
zm 

+ [ dr' ^(r',t)w{\r-r'\)i;(r',t)iP(r,t), (14) 



where the potential w(\r — r'\) describes the particle-particle interaction. We define the single-particle density operator 
via the usual definition, h(r,t) — xp' (r,t)ip(r,t) and the current density operator 

j(r,t) = -L [(V^{r,t))^{r,t)-^(r,t)V^{r,t)] . (15) 
2m 

It is lengthy but straightforward to show that these two operators follow the dynamics induced by the coupled 
equations of motion 

d t h{r,t) = -V-j(r,t) (16) 
md t ji(r,t) = ~h(r,t)ViV ext (r,t) -V jfi^(r,t) 

-V f M) / air' ^(r',t)X7 lW {\r-r'\)i/j(r',t)i/j(r,t) (17) 



1 The first term on the right hand side of Eq. 
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where we have defined the kinetic stress tensor operator 



fi,j{r,t) = ^- i [V l ^{r,t)V^{r,t)+V 3 ^{r,i)V l ^{r,i) 



(18) 



From the equations of motion for the operators, we get immediately the equations of motion for their expectation 
values 

d t n(r,t) = -V-j(r,t) (19) 
md t ji(r,t) = -n{r 1 t)\/ l V ext {r 1 t)~\I J {f l . J {r,t)) 

- fdr' p 2 (r,r',t)Viw(\r-r'\) (20) 



where p 2 (r,r',t) — (^(r, t)ip'(r', t)ip(r' , t)ip(r, t)}. Another rather lengthy and involved calculation allows us to 
write the force density due to the particle-particle interaction as a second-rank tensor, provided the interactions are 
negligibly small at the boundary of the integration volume in equation ([20]) . The result is [1J, l25| 



V y ir,, = -±V S I dy^-^>- I d\p 2 (r + Xy,r-(l-X)y,t) 



V fdv ViVj dw{ ^ f 1 
2 'J dV \y\ d\y\ J 

dr' p 2 (r,r',t)VM\r-r'\), (21) 

so that we arrive at the dynamical equation 

md t j(r, t) = -n(r, t)V,V ext (r, t) - V 3 P 4J (22) 

where we have defined 

Pi,j = W i , j + (f iJ ). (23) 

From here, by using the definition of convective derivative and re-scaling the particle momentum so that the stress 
tensor reads 

Pi,j = W l , J + (f itj ) -mnvi vj , (24) 

one obtains the equations of motion for the particle and current densities in a form identical to the single-particle 
equations of motion (p~2|> - (|T3j) . 

Like Eqs. (fl"2")) and (|13p . which, for a given initial condition, constitute a closedset, i.e., their solution is equivalent to 
the solution of the single-particle time-dependent Schrodinger equation, also their many-body counterpart, Eqs. (|19[) 
and (|22| are equivalent to the solution of the many-body time-dependent Schrodinger equation. This equivalence 
is a direct consequence of the theorems of time-dependent density- functional theory [53, [3] . These theorems state 
that, given an initial condition, there exists a one-to-one correspondence between the time evolution of the particle 
density n{r,t) and scalar potential V ex t(r, t) applied to the quantum mechanical system. A similar correspondence 



holds between the current density j(r,t) and an external vector potential A{r, t) [271 |29L | 30L l3l|. while the mapping 
does not generally exist between the current density and the external scalar potential [32j. The physical relevance 
of these theorems to our case is that the stress tensor Pij (r, t) in (|23|) is a functional of either the density or the 
current density, i.e., Pij(r,t) — Pij[n(r,t'),t] or Pij(r,t) — Pij[j(r,t'),t] (with t 1 < t). This implies that once the 
exact many-body stress tensor P^j is known, one can, in principle, recover from the solution of Eqs. (119|) and (|22p 
full information on the many-body wave-function. 

Needless to say, the exact stress tensor is unknown. However, starting from Eqs. (fT9|) and ([22|) one can develop 
perturbation schemes to approximate the exact stress tensor 3^, 33, 34|, at least for the problem at hand, thus simpli- 



fying enormously the solution of the many-body problem. In the following, we will describe one of such approximation 
schemes for the present case of current flow in a nanojunction. We will derive an equation of motion for the stress 

(3) 

tensor Pij and show that it depends on the so-called three-particle stress tensor P^j k , which in turn describes the 

(3) 

way three particles interact. The derivation of an equation of motion for P^j k would bring us into the maze of a 
hierarchic set of equations for stress tensors that describe electron-electron interactions to all orders. We will show, 
however, that for the case at hand, we can truncate this hierarchy and obtain a closed equation for the stress tensor 
P 
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VISCO-ELASTICITY OF THE ELECTRON LIQUID 

In parallel with the hydro-dynamic description of Quantum Mechanics, a visco-elastic formulation of the dynamics 
of the electron liquid has been derived within linear-response theory. It has been realized that a certain class of 
low-energy, long-wavelength excitations of the electron liquid may be mapped into the dynamics of a visco-elastic 
medium [351 ] . The dynamics of this visco-elastic medium is described by an equation of motion for the current density 
given by (in linear response and d dimensions, d > 1) 



mndtv(r,t) = 



f< I -(!-§)„ 



V (V • v(r, t)) 

f)V 2 v(r,t) - n{r,t)VV ext {r,t) (25) 



where K and fi are two complex constants which depend on the electron density n. These complex constants are 
expressed in terms of the more familiar viscosities, £ (bulk viscosity) and r\ (shear viscosity) and elastic constants K 
(bulk modulus) and [i (shear modulus) via the relations 

K{w)=K-iuQ, (26) 
= n — iuirj, (27) 

where ui is the frequency of the external perturbation used to excite the electron liquid. 

The next step is then to express the visco-elastic coefficients of the liquid in terms of its microscopic properties, 
i.e., relate these quantities to the resp onse functions. Here we only report the results that are relevant to the present 
work and refer the reader to Ref. [35j for an explicit derivation. We are only concerned with the DC (zero frequency) 
limit of the above quantities. By using an interpolation of the numerical results of mode-mode coupling theory [361 ] 
one finds the following density dependence of the zero-frequency shear viscosity (the bulk viscosity is identically zero 
in the same limit) [35| 



V 



n 60r s 3/2 + 80TS 1 - 40r s 2/3 + 62r7 1/3 

in 3D and in 2D by 



(28) 



r 2 2 N _1 

. — + 0.25r 2 



+21r7 2 + 23r 5 7 1 / 2 + 13, (29) 
where r s is the electron constant for the electron liquid with uniform density n: 



(3/47m) 1/3 3D 
(l/7m) 1/2 2D 



I Wit) OLJ , , 

r s a B = <i ,1/2 nn ' ( 30 ) 



with as the Bohr radius. 

It is interesting to point out that specific confining potentials (e.g., an electron liquid in a quantum well) may 
make the approximations used to derive Eqs. ([28]) and (129]) ill founded, leading to a peculiar behavior of the viscosity 
coefficients 1371. 



HYDRODYNAMIC APPROACH TO TRANSPORT IN NANOSCALE SYSTEMS 

In this section we show that in the case of nanoscale systems the stress tensor can be approximated to a form 
similar to the classical Navier- Stokes one. This is due to the geometric constriction experienced by electrons flowing 
in the nanostructure which gives rise to a very short "collisional" time [38L [3~ol ] . The system we have in mind is some 
nanoscopic junction sandwiched between two mesoscopic or macroscopic leads (see Figure [IJ and current is induced 
in the system by, e.g. polarizing the leads with a finite bias. In this regime, we show that one can truncate the infinite 
hierarchy of equations of motion for the electron stress tensor given in (|23|) to second order and thus derive quantum 
hydrodynamic equations. To realize how the simple presence of the junction has such a strong impact on the equation 



V 



FIG. 1: Schematic of a nano-junction connected to two bulk electrodes. A steady current is flowing from one electrode to the 
other. 



of motion of the current, one has to keep in mind that the former acts as a single impurity potential that cannot 
be avoided by the electron flow. This is different from the corresponding effect in bulk materials for which a certain 
density of impurities is necessary to have a finite resistance. 

Let us then employ the quantum Boltzmann equation for the single-particle distribution function f(r,p,t) (which 
can be derived from the time-dependent Schrodinger equation with standard techniques [4pj | ) and show how the short 
collisional time induced by the nanostructure allows us to close the equations for the stress tensor. 2 The quantum 
Boltzmann equation for the distribution function in a co-moving (Lagrangian) reference frame moving with velocity 



v{r,t) is [1J,|41| 



/[/] = D t f(r, P ,t) + *Lvf(r, P ,t) + eV^^M V 3%M 
m op op 

df(r,p,t) 

-mD t v (31) 

op 

where I is the usual collision integral pZj] , if is the sum of the external potential and the Hartree part of the interaction 
potential. The collision integral contains two terms, one elastic and the other inelastic. In what follows, it is important 
to realize that both terms can drive the system toward a local equilibrium configuration. 

From the quantum Boltzmann equation, we can derive the equation of motion for the moments of the distribution 
function. The general expression for the mth moment is the mth-rank tensor 

Ph,...,i m = I ' dpp h ...p im f(p,r,t). (32) 

The zeroth order is the single particle density, the first moment is the velocity field, and the second moment is the 
stress tensor we want to approximate. The equation of motion for the stress tensor contains a term proportional to 
the third moment P^: 

AP,j + PijV • v + P lM V kV] + Pk.jVkV, + VfeP^ = 



- / dp I[f]p iPj . (33) 
m J 



We note that P^ enters in (|33| only through its spatial derivative. If the latter is small then the hierarchy can 
be truncated [3, EJ . From ([33| we easily see that this derivative is small compared to the other terms whenever 
7 = ?i/(Lmax(aj, v c )) <C 1. Here u is the average electron velocity, L is the length of inhomogeneities of the liquid 
that give rise to scattering among three particles, u> the system proper frequency and v c the collision rate. The 
parameter 1/L enters through the spatial derivative of p( 3 \ lu from the frequency dependence of the interactions 
(in the DC limit of interest here u> — > 0), v c through the collisional integral /[/] oc —v c {J — /o), where /q is the 



2 Clearly, for the definition of local equilibrium distribution to be valid any length scale entering the problem has to be larger than the 
system Fermi wavelength. 
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equilibrium Fermi distribution. This derivative is indeed small for transport in nanostructures: When electrons move 
into a nano-junction they adapt to the given junction geometry at a fast rate, and approach to local equilibrium 
occurs at this fast rate even in the absence of electron interactions (38l . [3^ |. This "relaxation" mechanism occurs 
roughly at a rate v c — (At) -1 ~ (ft/ AE)^ 1 , where AE is the typical energy spacing of lateral modes in the junction. 
For a nano-junction of width I we have AE ~ ir 2 H 2 /m£ 2 and At ~ m£ 2 /Tr 2 h. If £ = 1 nm, v c is of the order 
of 10 15 Hz, i.e., orders of magnitude faster than typical electron-electron or electron-phonon scattering rates. The 
condition 7 = «/(Imax(u,i' c )) <C 1 thus requires the length of inhomogeneities L 1 nm, which is easily satisfied 
in nanostructures. Note instead that in mesoscopic structures this condition is not necessarily valid. In that case, the 
dominant relaxation rate v c is given by inelastic effects, i.e. it is of the order of THz, so that for typical lengths of 
mesoscopic systems, 7 ss 1 in the DC limit. Nonetheless, the above condition could still be valid for high-frequency 
excitations, like plasmons, and/or very low densities, so that moments of the distribution of order higher than two 
are negligible. 

By neglecting Vfci^J fc in (|33|) we can thus derive a form for P^-. Let us write quite generally the stress tensor 
Pi j as Pi j — SijP — TTij, where the diagonal part gives the pressure of the liquid, and m t j is a traceless tensor that 
describes the shear effect on the liquid. From (J33J) we thus find that the tensor ■n^j can be written as (in d dimensions, 
d>l) 



Ti,j = V ViVj + VjVi - -5 hj X7 k Vk ) (34) 



where 77 is a real coefficient (the viscosity) that is a functional of the density 4l| . We point out that ([34]) is in fact a 
particular case of a general stress tensor with memory effects taken into account [30I [3a . |42| . In our derivation this is 
the first non-trivial term of an expansion of the stress tensor in terms of the density and velocity field. Consequently 
the Navier-Stokes stress tensor in (I34| ca n be seen as the first-order (non-trivial) contribution to the exact stress 
tensor of the electron liquid (see also |25l. l30l. li^l). 

Using this stress tensor we finally get from (f22|) the generalized Navier-Stokes equations for the electron liquid in 
nanoscale systems 

D t n(r,t) = — n(r, t)V • v(r, t), 
mn{r,t)D t Vi{r,t) = -V t P{r,t) +V ^(r.t) (35) 
-n(r, t)ViV ext (r, t). 

Equations (|36|) are formally equivalent to their classical counterpart [l5j [see Eq. ([1])] and thus describe also nonlinear 
solutions, i.e., the possibility to develop turbulence in the electron liquid in its normal state. In the following, we 
will consider only the case in which the liquid is in the laminar regime and incompressible so that the viscoelastic 
coefficients are spatially uniform. This latter approximation is practically satisfied in metallic quantum point contacts 
(QPCs) but needs to be relaxed in the case of QPCs with organic/metallic interfaces (see, e.g., [2l|). In addition, 
for this case the Hartree potential is constant and its spatial derivative is thus zero. Therefore, (l36|) reduce to the 
Navier-Stokes equations for the density and velocity of a viscous but incompressible electron liquid 

D t n{r,t) = 0, 

V-u(r,t)=0, (36) 
mn{r, t)D tVl (r, t) = -V 4 P(r, t) + T]\7 2 Vl {r, t) - n(r, t)\7 z V ext {r, t). 



HEAT EQUATIONS FROM HYDRODYNAMICS 

The above results allow us to treat heat generation and transport using a simplified hydrodynamic approach. In fact 
we know that the flow of a viscous fluid, as described by our formalism, generates internal friction and consequently 
an effective temperature distribution inside the system. Therefore, when a steady state has been reached, we can 
supplement the Navier-Stokes equations with an equation for the energy balance. In the process of heat production, 
we need to identify a heat source, a mechanism for the dissipation of this heat and, since the system is in a steady 
state, equate these two terms with the local entropy production. In a recent paper [7[ we have developed this model 
obtaining the equation for the energy balance 



TTiji^djViir) + V ■ [A(r)VT e (r)] = c v (T e )v(r) ■ VT e (r), 



(37) 
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where T e is the electronic temperature, k(r) is the diffusion constant and cy is the specific heat at fixed volume of 
the electron gas. 3 Eq. ([37)) can be either justified on physical grounds, or derived formally as high-order expansion 
of the many-particle stress tensor [25| . We also stress once more that in deriving this equation we have assumed that 
the flow of the electron liquid is laminar, i.e., we are far from the onset of a turbulent regime 14, 15j]. Obviously, in 




writing Eq. (|37[) we have assumed that some thermodynamic quantities like temperature and entropy for an electron 
liquid flowing in a nanostructure can be defined. This is a much debated point, and obviously we do not have a general 
solution for it. However, here we argue that the electron temperature may be defined as the one ideally measured by 
a probe weakly coupled to the system and in local equilibrium with the latter While this operational definition 
may not be simple to realize in practice, we know from experiments that local heat generation due to current has a 
large effect on the stability of nanostructures [9( . 

From the form of Eq. (I3T[) we can deduce a general relation between the applied bias and the electron temperature. 
To do this, we realize [7\ that the electron fluid velocity, v, (which is generally smaller than the Fermi velocity [22j]) 
responsible for the transport of current and heat is, in linear response, proportional to the bias V . 4 This simple 
proportionality, and the usual result that k oc cy, bring us to the general result 

T e = J ee V, (38) 

where j ee is a constant whose expression in terms of microscopic parameters of the electron liquid has been recently 
derived for a quasi-adiabatic connection between the leads and the nanojunction Q 

^ = L16x (-^)\l z ^rz (39) 

where G is the conductance of the system in the limit of zero bias, A c its cross section, d is the dimensionality (d > 1). 
Moreover, 

7 = k 2 F k 2 B X e /9 (40) 

in 3D, and 7 = irkpk^Xe/G in 2D [7j, kp is the Fermi momentum, ks the Boltzmann constant, and A e is the inelastic 
mean free path. 

Interestingly, Eq. (|38[) can be obtained from general thermodynamic arguments, by comparing the energy dissipated 
in the transport process in the nanostructure (proportional to V 2 from Ohm's law), and the energy carried away by 
electrons (proportional to T 2 for small temperatures) [7j . 

LOCAL ELECTRON HEATING 

In the case of a finite background temperature and in the absence of ionic heating, from our hydrodynamic theory, 
the local temperature of the electrons in the nanostructure is given by [3] 

T e (V) = sfn^y 2 (41) 

where V is the external bias, and To is the electron temperature deep into the electrodes. If we now let the ions heat 
up, their effective local temperature is given by 0| (for values of the parameters such that the argument in the root 
is non-negative) 

T={Tt + 1 i p V 2 -ltV i ) l/i (42) 

where 7 ep can be expressed in terms of the physical parameters of the nanostructure [H, l43| , and we have assumed 
that both the ions and the electrons are at the same background temperature To deep into the electrodes. At zero 
background temperature and for negligible electron-electron interactions from the above equation we obtain the known 
result for the local ionic temperature [3|, 0, 51] 

T ~ j ep W. (43) 



3 For an electron gas at low temperature, cy = cp since the correction is second order in temperature. 

4 V may be given by an external battery, or the potential due to a charge imbalance. 
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Effect on conductance - We can now calculate the effect of local electron heating on the conductance of a nanos- 
tructure. We focus on the quasi-ballistic regime and we generalize Eq. (13) of Ref. [43| for the inelastic current in 
the presence of a finite electron temperature. 5 We also consider one mode frequency w. We will generalize later to 
more modes. To take into account the effect of an effective local electron temperature on the inelastic current, one 
faces the calculation of terms with factors of the type dE f%{\ — fgj-J), with a, /3 = {R, L} corresponding to 
electrons moving from either left or right, and /g = (exp((E — ^ Q )/fcBT) + l)^ 1 is the Fermi distribution with the 
difference between the electrochemical potentials equal to the bias, fJ>L — (j,r = V. (Refer to [43| for additional details 
on the notation.) 

We could provide a numerical calculation of the inelastic current. However, we are interested in an analytical 
expression and thus proceed as follows. We evaluate the above integrals in the Sommerfeld approximation and keep 
only the terms of zeroth order in the electron temperature (this is reasonable because the local electron temperature 
is generally a small quantity). This approximation brings us to the expression for the current flowing in the system 



I : 



i/4 log(e^ v -") + 1) 



(44) 



where 7/ is the amplitude of the conductance drop at V = V c = u) for zero electron and phonon temperature, G e i is 
the elastic conductance at zero bias, and (3(V) — l/ksT e (V) where fcs is the Boltzmann constant. By differentiating 
Eq. (|44p with respect to bias, and again keeping only the terms of zeroth order with respect the electron temperature, 
we arrive at 



G k B 
7J- - 1 li 

G e l V 



T, 



1/4 



1 



2 /3(V)(w-V) + 1 



(45) 



To obtain this result, one also has to bear in mind that the approximations we make pertain to the energy region 
V ~ ui, thus d(/3(V)(V — u>)/dV ~ 0(V)- An expression for the conductance similar to Eq. (|45p can be derived for the 
case of zero electron temperature [43j, i.e., (3 — ► 00. Notice, however, that for consistency, one has to take this limit 
in the expression for the current (l44|) before taking the derivative with respect to the bias. 




FIG. 2: Plot of G/Go as given by Eq. (HJ where Go = 2e 2 /h. The solid (red) curve has been generated with the following 
parameters: j ee = 180 K/V, 7 ep = 600 K/VV, G e i = 0.98 x Go, To = 5 K ,71 = 0.02, hcu = 0.05 eV. The dashed (green) line 
corresponds to the case of zero phonon and electron temperature (j ee = y ep = 0). The dashed-point (black) curve corresponds 
to the case when only the phonon heating is taken into account (7 ee = 0). In the figure we also define the parameter A (see 
text). 

An example of the effect of local ionic and electron heating on conductance is given in Fig. [2] (see also discussion 
below) . In the absence of both effects (and at zero nominal background temperature) the conductance shows a simple 
step-like drop at the bias corresponding to the energy of the phonon mode. The ionic heating introduces a shoulder 



5 Note that a factor 2 is missing in Eqs. (7) and (8) of Ref. fijjll . 
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at biases larger than the mode energy, while the electron heating broadens the conductance curve with an effective 
temperature larger than the nominal background temperature. 

Comparison with experiments - To compare our results with available experimental data we consider a D 2 molecule 
sandwiched between two Pt leads [l3| . We focus on the predictions of our hydrodynamic theory on the local electron 
heating effect. Therefore, we do not attempt to do a full first-principles calculation of ionic heating, and take the 
relevant parameters from experiment. For the D 2 molecule we consider a cross section of tt x l A 2 , i.e., a circle 
with radius ~ 1 A. The nominal electron temperature deep inside the electrodes is taken to be To = 5 K. From the 
experimental results we have the frequency of the phonon mode V c — 0.05 eV, the drop of the conductance 



and the conductance at zero bias, G e i = 0.984 Go- We use as fit parameters 7 ep and evaluate 7 ee from Eq. (|39|) . In 
obtaining j ee we have assumed an inelastic mean free path A e of 1 /im, a value in line with the expectations for this 
system [Hj]. We have also assumed that the electron density that enters the local heating is the one of the chemical 
bonds between the D and Pt atoms. This density is estimated to be close to the Pt bulk density, n = 6.6 x 10 28 m~ 3 
which gives the electron constant r s ~ 3. From these values, the electron viscosity r/ and the constant j ee are easily 
obtained from Eqs. and (|4"0")) . respectively: The electronic heating constant is predicted from Eq. (f3"9"|) to be 
7 ee = 180 K/V. This implies an effective electron temperature of about 10 K at the D 2 junction at a bias of 50 mV. 
This temperature is higher that the nominal bulk temperature. 

The ionic heating constant is found to be 7 ep = 405 K/vF. 6 This value can be compared with the corresponding 
7e P for a Au point contact at small biases which is about 7^" = 170 K/y/V [43j. This means that the Pt-D 2 - 
Pt system heats up more than the Au QPC. For instance, at 0.1 V the ions of the Pt-D 2 -Pt junction have an 
average temperature of about 130 K while at the same bias the gold atoms heat up locally to about 54 K. This 
larger temperature is reasonable since, while the conductance is similar for a Au point contact and Pt-D 2 -Pt, the D 2 
molecule is lighter than Au with a consequent increase of the electron-phonon coupling. In addition, the modes of 
the D 2 molecule have lower probability to elastically scatter into the bulk modes of Pt - thus reducing lattice heat 
dissipation into the bulk electrodes - than the modes of a single Au atom into the bulk modes of Au. Both effects 
lead to a higher local ionic temperature. We thus expect the Pt-D 2 -Pt junction to be more unstable under the same 
bias conditions than a Au point contact, i.e., we expect that the chain Pt-D 2 -Pt breaks, on average, at much smaller 
biases than Au point contacts due to heating effects. 

The theoretical conductance containing both the local electron and ionic heating effects is reported in Fig. [31 
together with the experimental data. The qualitative agreement between theory and experiment is very good. It 
is interesting to note that the tail of the experimental data goes approximately as V ' 8 , while the theory predicts 
It is important to realize, however, that at large biases, other effects such as current-induced 
forces and other structural instabilities may also contribute to the actual value of the conductance [HI]. 



li = 



d(I-Ig) 
dV 



= 0.02, 



(46) 



jig 
dV 




0.98 



o 



o 

£5 0.97 



0.96 







0.02 



0.04 



0.06 



0.08 



V[V] 



6 This confirms that the ions heat more than the electrons at the same bias and our approximation that leads to Eq. (142 [ 1 is justified. 
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FIG. 4: A as a function of vibrational mode energy as defined in the text. We have tested the linear behavior for different 
values of the percentage a at which we calculate A. The lines that connect the symbols are linear regressions. 



FIG. 3: Comparison between the experimental data [TJ and our theory [solid line, (red)] [l4[. 

We have also performed a second fit, not shown here, using 7 ee and 7 ep as free parameters. The values for these 
parameters obtained from this second fit are close to those obtained from the theory and the one-parameter fit by 
less than 10 % (we find the best fit for 7ee = 200 K/V). 

Inelastic conductance width - Let us now discuss how the width of the inelastic conductance around the vibrational 
mode increases with bias (see Fig. [2j. This quantity can be directly measured and provides additional information 
on local electron heating. If the background temperature is zero, the local electron temperature increases linearly 
with bias as in Eq. (j41"j) . Let us define the quantity A as shown in Fig. [2j It is the energy distance between the 
middle drop of the conductance and the value at which the conductance assumes its purely elastic value within a ratio 
a = [G e i — G(V c )]/G e i as indicated in Fig. O This quantity is plotted in Fig. [1] for different values of the vibrational 
mode energy and for a few values of a, assuming that the vibrational energy is the only quantity allowed to vary. We 
conclude that the width A increases linearly with the vibrational energy to reflect the linear bias dependence of the 
local electronic temperature. A systematic experimental study of this quantity would thus provide more information 
on the electron heating phenomenon. 



DISCUSSION 



Our analysis in conjunction with the experimental data suggests that electrons heat up locally at the Pt — D2 — Pt 
junction. Our Eq. (|41[) also predicts that electrons cool down when lowering the bias. On the other hand, a constant 
electron temperature - above the background temperature - for all biases is difficult to understand on physical 
grounds, unless one assumes the existence of an external source of energy that keeps the electron hot even at zero 
bias. 

Experimental data showing an electron temperature equal to the background temperature, i.e., negligible electron 
heating, may be consistent with the fact that the effective cross section "seen" by the electron liquid is the one of a Pt 
atom and not of a D2 molecule. 7 If that were the case, the effective cross section would be 7 times larger than that 
of the D2 molecule, and since the electron temperature scales inversely proportional to the cross section (see Eq. (|39|) 
and Ref. [14j ) the electron heating temperature would be lower than the background temperature. The conductance 
on the other hand is unlikely to be so sensitive to this cross section due to the extended nature of the Pt d-orbitals. 

Further generalization of Eq. (|45[) to the case where many vibrational modes are present is possible. For example, it 
has been reported that a H2 molecule sandwiched between two Pt leads shows two fundamental vibrational frequencies 



at 48 meV and 62 meV |13|. If one assumes that scattering by these two modes is uncorrelated, a straightforward 



7 This system can be thought of as a junction Pt-Pt-Pt with just one or few platinum atoms forming an effective QPC which is not 
significantly affected by the presence of the deuterium molecule. 



12 




i i i i i i i i i i 

0.02 0.04 0.06 0.08 

V[V] 

FIG. 5: a) Experimental results of the non-linear DC conductance of a H2 molecule between two electrodes. The two steps 
coming from the two phonon modes at 48 and 62 meV are clearly visible, together with some unexpected structure at low bias, 
b) Plot of Gh 2 as a function of the external bias according to Eg. 1471 In this plot, that is not a fit to the experimental data, we 
have used 711 = 0.01, 7/2 = 0.002, 7 ee = 180 K/V, 7 ep ,i = 7e P ,2 = 400 K/Vv. The other parameters are the same as in Fig. [3] 



generalization of Eq. (|45[) leads to 

Gh 2 ^ 

Gel 



1 7/1 {T +7ep,l^ 



7/2 (T + 7 2 

U>2 



V 2 



-Tin V4 
7i^ 4 ) 1/4 



1 



e 0(V)(^-V) + 

1 

z [i{V)(u 2 -V) + I 



(47) 



where wi and u>2 are the two vibrational frequencies and we have taken into account the possibility that the two 
coupling constants 7ep,i and 7ep,2, and the two amplitudes of the conductance drops 771 and 772 be different. A plot 
of Gt7 2 is reported in Fig. [5] as a function of the external bias along with the experimental data. Since the cross section 
for D2 and H2 is essentially the same, and the electron heating does not depend on the mass of the ions, our estimate 
of 7 ee holds for H 2 as well. Our results are again in qualitative agreement with the available experimental data [l3j ]. 
although our theory might be not sufficient to quantitatively describe all the experimental findings. Indeed, our fit 
in this case has failed in producing any sensible result for the constants 7 ee , 7 ep ,i and 7 ep , 2 : the large fluctuations of 
the experimental data, especially in the region of small bias and close to the phonon modes energies do not allow for 
a systematic fit of the data with the theory. Finally, it is interesting to note that a value of 7 ep similar to the one we 
have obtained for the D2 molecule gives a reasonably good agreement between theory and experiment also for the H2 
molecule. This seems to suggest that the longitudinal modes of the bonds between the H and Pt atoms are mainly 
responsible for the local ionic heating of the Pt-H2-Pt junction, and similarly the longitudinal modes of the bonds 
between the D and Pt for the Pt-D2-Pt junction. We expect that such modes are slightly affected by the change of 
mass of the smaller atom in the bond. Clearly, more theoretical and experimental work in this direction is necessary. 



CONCLUSIONS 



We have discussed a novel hydrodynamic approach to transport that allows the description of charge and heat 



flow in terms of the single-particle density and velocity field of the electron liquid [14j . The theory allows us to make 
predictions about the electron flow past a nanostructure and its dependence on the external bias (or the current). One 
such prediction is the heating of electrons locally at the nano-junction [7J. Here we have considered the measurable 
consequences of this effect on the inelastic conductance which shows a broadening at the inelastic step larger than the 
one expected from the background nominal temperature. We have compared our theory with available experimental 
results [lH and found a reasonable quantitative agreement for the case of a D2 molecule between two Pt leads. For 
the case of a H2 molecule between the same leads our theory is only in qualitative agreement with the experimental 
findings. We also predict that the width of the inelastic conductance step should increase linearly with bias, a fact 
that can be tested experimentally. 
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